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1. Introduction 

There are two well-known problems which are related to each other. They are the 
discrepancies in gravitational constant experimental data and the hierarchy problem. 
Discrepancies (see e.g. Figure 2 in the "CODATA Recommended Values of the 
Fundamental Constants: 2006") are usually explained by extreme weakness of gravity. 
It is very difficult to measure the Newton's gravitational constant Gn- Certainly, for 
this reason geometry of an experimental setup can effect on data. However, it may 
well be that, the discrepancies can also be explained (at least partly) by underlying 
fundamental theory. Formulas for an effective gravitational constant following from 
such theory can be sensitive to the geometry of experiments. For example, if correction 
to the Newton's gravitational potential has the form of Yukawa potential, then the 
force due to this potential at a given minimum separation per unit test-body mass 
is least for two spheres and greatest for two planes (see e.g. [1]). Therefore, an 
effective gravitational constant obtained from these formulas acquires different values 
for different experimental setup. 

The hierarchy problem - the huge gap between the electroweak scale Mew ~ 
lO^GeV and the Planck scale Mp;(4) = 1.2 x lO^^GeV - can be also reformulated in 
the following manner: why is gravity so weak? The smallness of Gn is the result 
of relation Gn — Mp^^^^ and huge value of A/p;(4). The natural explanation was 

proposed in [21 13] : the gravity is strong: Gv — ^^pi^-i))'' ^ ^^ew^'^^ ^^'^ happens 
in (T) = 4 + d)-dimensional spacetime. It becomes weak when gravity is "smeared" 
over large extra dimensions: Gjv ~ G-p/Vd where Vd is a volume of internal space. 

To shed light on both of these problems from new standpoint we intend to 
investigate multidimensional gravity in non-relativistic limit. To do that, first, we are 
going to obtain a solution of (Z? = 3 + d) -dimensional Poisson equation in the case of 
toroidal extra dimensions. From school years we know that newtonian gravitational 
potential of a body with mass m has the form of ^plr^) — —GNtn/r^ where Gn is 
Newton's gravitational constant and — Ir^l is magnitude of a radius vector in 
three-dimensional space. This expression is the solution of three-dimensional Poisson 
equation with a point-mass source (and corresponding boundary condition fir^) — > 
for 7'3 +oo) or it can be derived from Gauss's flux theorem in three-dimensional 
space. To investigate effects of extra dimensions, it is necessary to generalize the 
Newton's formula to the case of extra dimensions. Clearly, the result depends on 
topology of investigated models. We consider models where D-dimensional spatial part 
of factorizable geometry is defined on a product manifold Md = R'^ x T'^. M.^ describes 
three-dimensional flat external (our) space and T'' is a torus which corresponds to d- 
dimensional internal space. Let Vd be a volume of the internal space and a ~ V]^^"^ 
is a characteristic size of extra dimensions. Then, Gauss's flux theorem leads to the 
following asymptotes for gravitational potential (see e.g. [3j and our Appendix for 
alternative derivation): ip ^ l/ra for >> a and ip ^ l/rg^J^ for r^+d << a where 
and rs+d are magnitudes of radius vectors in three-dimensional and (3-|-(i)-dimensional 
spaces, respectively. Obviously, an exact solution of _D-dimensional Poisson equation 
should show smooth transition between both of these asymptotes. This formula 
gives possibility to investigate characteristic features of multidimensional gravity in 
non-relativistic limit which enable us to reveal extra dimensions (or to establish 
experimental limitations on extra dimensions). 

In our paper we obtain such exact solution for arbitrary number of the extra 
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dimensions. In the case of one extra dimension, this expression acquires very compact 
and elegant form. Found exact solution is applied to a number of practical problems, 
e.g. to calculate gravitational force between two spherical shells (or balls). It 
gives a possibility to calculate an effective four-dimensional gravitational constant 
for given configurations. For example, in the case of two balls we show that, in 
principle, the inverse square law experiments enable us to detect a deviation from 
the Newton's gravitational constant. Then, we generalize our model to the case of 
smeared extra dimensions. It means that we suppose that test bodies are uniformly 
smeared/spreaded over all or part of extra dimensions. We prove that the gravitational 
potential does not feel the smeared extra dimensions. In particular, if all extra 
dimensions are smeared, then the inverse square law experiment does not show any 
deviation from the ordinary Newton's formula, and this conclusion does not depend 
on sizes of the smeared extra dimensions. Nevertheless, the hierarchy problem can be 
solved in these models. 

The paper is structured as follows. In section 2 we get the exact solution of 
multidimensional Poisson equation in the case of spacial topology x T'^. This 
formula is applied to some practical problems in section 3 to get gravitational 
potential, gravitational acceleration of a point mass and gravitational self-energy for 
these problems. In section 4 we investigate gravitational interaction of two spherical 
shells. Then, in section 5 we generalize our model to the case of smeared extra 
dimensions. Here, we prove that gravitational potential does not "feel" the smeared 
extra dimensions and demonstrate that the hierarchy problem can be solved in this 
case. A brief discussion of the obtained results is presented in the concluding section 
6. In Appendix A, we get the expressions for the gravitational force law at small 
and big separation between point masses on multidimensional manifold with topology 
R3 X T"*. 



2. Multidimensional gravitational potentials 

In Z3-dimensional space, the Poisson equation reads 

Ad^pd = SoGvPoivD) , (2.1) 

where So = 2TT°^^/r{D/2) is a total solid angle (square of (£> — l)-dimensional sphere 
of a unit radius), G-p is a gravitational constant in [V ^ £> + l)-dimensional spacetime 
and pd{yd) = m5{xi)5{x2)...5{xD)- 

2.1. Spatial topology 

In the case of topology MP , (j2.1|) has the following solution: 

(^~2)7- 

This is the unique solution of (|2.ip which satisfies the boundary condition: 
lim fDi^o) — 0. Gravitational constant in (j2.ip is normalized in such a 

way that the strength of gravitational field (acceleration of a test body) takes the 
form: —dipu/dru = —G-pm/r^^^ . 



^d{vd) = - ,^ ' ^>3. (2.2) 
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2.2. Spatial topology M.^ x T"^, d-arbitrary 



If topology of space is M'^ x T'^, then it is natural to impose periodic boundary 
conditions in the directions of the extra dimensions: (/s^ (ra , , ^2 , ■ ■ ■ ,£,i, - ■ ■ ,£,d) = 
<PD(r3, Ci, ^2, ■ • ■ J + ctj! ■ • ■ ; Cd), i — where Oi denotes a period in the 

direction of the extra dimension ^i. Then, Poisson equation has solution (cf. also 
with [3l|4]): 

Givm 



Voir 3, £.1 



+ 00 
/C]^ ^ — c>o 

/27rfci 



+ 00 

E 

k^ — — oo 



X COS 



V 01 



6 



exp 



. cos 




(2.3) 



cos ( ^^Ci ) s-nd, second, 



To get this result we, first, use the formula (5(^i) = ^ J2k^- 
put the following relation between gravitational constants in four- and 2?-dimensional 
spacetimes: 

Sd Gt> 

0. In this limit zero modes 
ki = give the main contribution and we obtain v^nfr^, £i , — > —GMfn/fa- 



(2.4) 



U,=i at 

The letter relation provides correct limit when all 

give the main contribution and we obtain (Pd{'<^3^ ■ 
was widely used in the concept of large extra dimensions which gives possibility to 
solve the hierarchy problem [2l[3] . It is also convenient to rewrite (|2.4p via fundamental 
Planck scales: 



Sd 



M 



P((4) 



= M 



2+d 
Pl(V) 



(2.5) 



where MpK^^^ — G^^'^ — 1.2 x lO^^GeV and Mpi(^x>) = G^^'^""^"^ are fundamental 
Planck scales in four and 2? spacetime dimensions, respectively. 

In opposite limit when all —>■ +oo the sums in (12. 3p can be replaced by 
integrals. Using the standard integrals (e.g. from [5]) and relation (|2.4p . we can 
easily show that, for example, in particular cases d = 1,2 we get desire result: 
V.i,(r3, Ci, . . . , U) - ~Gvm/[{D - 2) rl+j]. 

Prom (|2.3p . it follows that potential energy of gravitational interaction between 
two point masses m^"''^ and to^''^ with radius vectors r^'' and r^' reads 



'l/(2+d) 




X COS 



(2.6) 
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2.3. Spatial topology x 

In the case of one extra dimension d — 1 we can perform summation of series in 
To do it, we can apply the Abel-Plana formula or simply use the tables of series [S]. 
As a result, we arrive at compact and nice expression: 

sinhf22La) 



¥'4(r3,C) = 



cosh (-n:^)- cos 



(2.7) 



where G [0, +00) and ^ e [0, a]. It is not difficult to verify that this formula has 
correct asymptotes when » a and r4 << a. Figure 1 demonstrates the shape of 
this potential. Dimensionless variables ryi = 7'3/a G [0, -|-cxd) and 772 = £,/a G [0,1]. 
With respect to variable 772, this potential has two minima at 772 =0,1 and one 
maximum at 772 = 1/2. We continue the graph to negative values of 772 G [—1,1] to 
show in more detail the form of minimum at 7/2 = 0. The potential (j2.7p is finite 
for any value of r3 if ^ 7^ 0, a and goes to —00 as — I/t'I if simultaneously ra — > 
and ^ — > 0, a (see Figure 2). We would like to mention that in particular case ^ = 
formula ()2.7p was also found in [6]. 




Figure 1. Graph of function ip(r]i,r]2) 

— sinh(27r?;i)/[7;i(cosh(27r»)i) — cos(27r7;2))]. 



'/'4(r3,C)/(G'jvni/a) 



2.4. Yukawa approximation 



Having at hand formulas (12. 3p and (j2.7p . we can apply it for calculation of some 
elementary physical problems and compare obtained results with known newtonian 
expressions. For a working approximation, it is usually sufficient to summarize in 
(|2.3p up to the first Kaluza-Klein modes |fci| = 1 (i = 1, . . . , c?): 



1 



2 ^ exp 



27r 



ra cos 



2tt 



-6 



(2.8) 



Then, the terms with the biggest periods a.; give the main contributions. If all test 
bodies are on the same brane (^^ ~ 0) we obtain: 



tpoirs,^! = 0,^2 = 0, = 0) = ipoira) 



r3 



1 + a exp 



(2.9) 
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Figure 2. Section ^ = of potential 112.71 1. Solid line is (^(r;i,0) = 
— sinh(27r7yi)/[»yi(cosli(27rr)i) — 1)] which goes to —l/rji (dotted line) for rji +oo 
and to —l/{nrif) (dashed line) for r]i — » 0. 



where a — 2s (1 < s < d) , X — a/(27r) and s is a number of extra dimensions 
with periods of tori a.; which are equal (or approximately equal) to a = max a.;. If 
ai = a2 = ■ ■ ■ — ad = then s = d. Thus, the correction to Newton's potential has 
the form of Yukawa potential. It is now customary to interpret tests of gravitational 
inverse-square law (ISL) as setting limits on additional Yukawa contribution. The 
overall diagram of the experimental constraints can be found in [7] (see Figure 6) and 
we shall use these data for limitation a for given a. 

3. Application {^i — 0) 

Now, we apply formulas (|2.3p and (|2.7[) to some particular geometrical configurations. 
For our calculations we shall use the case oi = £,2 — ■ ■ ■ — £,d = 0- It means that 
test bodies have the same coordinates in the extra dimension. It takes place e.g. when 
test bodies are on the same brane. Also, to get numerical results we should define the 
sizes tti of the extra dimensions. The e'^e~ leptonic interaction experiments at high- 
energy colliders show that there is no deviations from Coulomb's law for separations 
down to lO^^^cm [Tj. Therefore, if the standard model fields are not localized on the 
brane, this value can be used for upper bound of : ^ 10~"'^^cm (in section 5 we 
demonstrate how to avoid this argument). These values can be greatly increased 
if we suppose that the standard model fields are localized on the brane. In this case, 
we can obtain the upper bound for a — maxa^ from the gravitational inverse-square 
law experimental data depicted in Figure 6 of paper 7 . For example, the Yukawa 
approximation (j2.9p shows that a ^ 2 for d = I. For this value a. Figure 6 gives 
A < 4.7 X lO-^cm ^ a < 3.0 x lO^^^m. 

3.1. Infinitesimally thin shell 

Let us consider first an infinitesimally thin shell of mass m = AnR^a, where R and a 
are radius and surface density of the shell. Then, gravitational potential of this shell 
in a point with radius vector (from the center of the shell) is 



rs+R 




(3.1) 
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r-3 



k 



Xk 



sinh(27rxfc-R) , > R 



and 



R 



R+ra 
2nGM(TR f 

J Z^exp 

R-ra ^ 
1 



2Gn(jR 



~2TTXkr']dr' 



(3.2) 



^■3 



fe 



where 



Xk 



E 



1/2 



E 

A; 



+ 00 

E 

ki— — o 



-\-oo 

■ E 

kd = -o 



(3.3) 



and the prime in sums denotes that the zero mode fci = . . . = fc^ = is absent in 
summation. It means that the Newton's expressions (which correspond to the zero 
mode) are singled out from sums. In the case d = 1, these expressions can be written 
in the compact form: 



, , Gno' Ra 

</'4(?'3) = hi 



cosh 



^ 27r(r3+fl.) ^ 



1 



r3 



rs^R- 



(3.4) 



cosh(-M)„l 

These formulas demonstrate two features of the considered models. Firstly, we see 
that inside (ra < R) of the shell gravitational potential is not a constant. Thus, a test 
body undergoes an acceleration in contrast to the newtonian case (and to Birkhoff's 
theorem of general relativity in four-dimensional spacetime which states that the 
metric inside an empty spherical cavity in the center of a spherically symmetric system 
is the Minkowski metric). Secondly, if — > R, these potentials have a logarithmic 
divergency of the type: J^t^i ^/k. For example, in this limit R (|3.4|) has the 

following asymptotic behavior: 



<P4(?'3) 



Gntti 
R~ 



1 



2ttR 



■In 



27r|i?-r3| 



Gnttl 
R 



[1 + 5] 



(3.5) 



where we took into account R >> a and — raj << a. In particular case 2t:R — 10cm 



and 2'k\R — rj\ = 10 a, the deviation (5 constitutes 2.3 x 10 



"3 and 2.3 x 10"^° parts 
of the newtonian value —GM'm/R for a = 10~^cm and a = 10~^^cm, respectively. In 
principle, the former estimate is not very small. However, it is very difficult to set 
an experiment which satisfies the condition \R ^ ry,\ << a. If the shell has a finite 
thickness, then the divergence disappears. 



3.2. Spherical shell 

Here, we consider a spherical shell of inner radius Ri and outer radius i?2 and mass 
m = 47r p{R)R^dR where p{R) is a volume density of the shell (in the case of 

constant volume density p = m/ [^(i?2 " -^i)])- Then, the potential outside of the 
shell is 

{-R.2 
y]_^-2.x.r, f p(^R)Rsinh{27TXkR)dR 



(3.6) 
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Gntti GmP 



1 



27r2r3 



^ — exp[-27rxfer3]%(i?)|^J , (3.7) 



Xk 



where 

hk{R) = 2TTXkRcosh[2TTXkR] - sinh[27rx/c-R] ■ 
Inside of the shell we obtain: 



(3.8) 



<PD{r3) = — sinh(27rxfer3) / p(i?)i?e-2-««di?^ 

''3 . Xfc J 



and for the constant p: 

foirz) = -^ttGnP {Rl - Rl) 



1 



-2-!TXkR 



R2 

Ri 



(3.9) 



(3.10) 



To get potential within the shell, we can use the following relation: 

^D(r3) = (/p^'Hr3|i?2 = r3) + <pgHr3|i?i=r3), i?i < < i?2. (3.11) 
Thus, in the case of constant p the gravitational potential within the shell reads 



Voir 3) = 2770 nP 



i?2 



2^\ 
3r3 y 



GnP v^' J_ 

2^ ,.3 



27r2r3 



X [27rxfer3 - sinh(2^Xfe^3) (2^Xfci?2 + 1) e-2-x.fl2 _ e-^^^-^-^'hkiRi)] . (3.12) 
For one extra dimension d = 1 we obtain a compact expression which is valid for full 
range of variable rs > 0: 



R2 ( I, f 2Tr(R+r 

GNpa /• J cosh(^^^ 



^3 



/ 



i?ln< 



- 1 



lcosh(^ 



dR, 



1 



(3.13) 



where p is taken to be constant. 

It can be easily seen that all these potentials arc finite in the limit — > i?i , R2 
and Ri ^ i?2, i-c. divergency of the potential is absent for finite thickness of the shell. 
However, divergency takes place for acceleration of a test body. We can see it from 
exact formulas for acceleration outside of the shell: 



drs 



GnP 
2-k'^tI 



^ (27rxfcr3 + 1) 6-2-^'='-=* {hk{Ri) - hk{Ri)] < , rs > R2, (3.14) 



within the shell: 



— — = - ttGnP — 
drs 3 



RI 



GnP 
2ttHI 



Y! \ {hkira) {27rxkR2 + 1) e-^^^'""- 



Xk 



- (27rxfer3 + 1) e-^^'^'^-^hkiRi)} , Ri<r3<R2 



(3.15) 
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and inside of the shell: 

dipD GnP v^' 1 , , n 
T — "fc(^3 



R 



27rxfe 



^0, r3<i?i. (3.16) 



All of these equations p.l4p - p.l6p have logarithmic divergency in the limit 
i?2, e.g. in d = 1 case: 



difi 
drs 



— — — > ^2GNpa\n 



27T\Ri.2~r3\ 



(3.17) 



where sign "-" corresponds to Ri and sign "+" corresponds to — > i?2- 

Putting the limit Ri = 0, i?2 = i? in expressions obtained above in this subsection, 
we can get the corresponding equations for a sphere. For example, using the same 
conditions as for (|3.5p . we can get for a sphere {d — 1): 



dip4 
drs 



1 



3a , f2TT\R- 

In 

2'kR \ a 



R 



[1 + 6] 



(3.18) 



Therefore, deviation from the newtonian acceleration S — 36 and we can conclude that 
this deviation is also difficult to observe at experiments for considered parameters. 

(|3.14p and (|3.16p show that acceleration changes the sign from negative outside 
of the shell to positive inside of the shell (see Figure [3|). This change happens within 
the shell. 
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Figure 3. Graph of acceleration w = —dipi/drs in dimensionless units (see 
H3.14|l - H3.16|l '). Here, a = 2.5, Ri = 5 and R2 = 10. The dashed lines correspond 
to radii of the shell. The rightmost line goes to the newtonian asymptote — 1/rg 
when r3 — > -|-oo. Point ra = corresponds to unstable equilibrium. 



3.3. Gravitational self-energy 



Gravitational self-energy of the spherical shell with potential lpd{tz) given by formulas 
(|XT^ and ((XT^ reads 



U = 2Tip I r;^ipD{r3)dr3 

jRi 

QGntu^ 



SGMm^ 



( — 2i?2 — 3i?J + dR^R^^ 



10 (i?i - Rf) 
;{Rl~Rl) + 



:Ek 



(27rxfe)' 



(3.19) 
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where 

Ek = il + 27rxfci?2)e-'"^'=^^ [-hk{R2) + 2/ifc(i?i)] 

- (1 + 2TTXkRi)e-^'''"''^'' hk{Ri) . (3.20) 

In the case of sphere i?i = 0, i?2 = -R and this expression is reduced to the 
following formula: 

^^_3G^r 15 



X V' — 

k 



5R { 2{2ttR) 

2TTXkR f 1 



3 



2 

{l + 2TiXkR)e-^''''''''hk{R) 



3 \2^XkR, 

= Un + 6U. (3.21) 

Here, Un and 6U are the newtonian self-energy of the sphere and deviation from it, 
respectively. In the case of one extra-dimension d = 1 and for condition a ^ R we 
obtain: 

5 / a \2 1 5 / a \2 ^ ^ 

^^-^-4^y Ep-^-^y ^ (3-22) 

where we used that J^t^i ^/k'^ — t^'^ It is worth of noting that the correction to 
the newtonian formula has power law suppression instead of exponential one as it is 
usually expected in Kaluza-Klein models. Nevertheless, for astrophysical objects this 
correction term is negligibly small because of big difference between a and R. For 
example, if a = lO'^cm, we get that 5[/ « 4 x IQ~'^''Un and 5U k. Ax IQ'^'^Un for 
Sun (i? w 7 X 10^°cm) and neutron star {R 7 x lO^cm), respectively. 

4. Gravitational interaction of two spherical shells 

Let us consider two spherical shell with radiuses Ri < R2 and constant density p 
for first shell and radiuses R[ < R'2 and constant density p' for second shell. Then, 
potential energy of gravitational interaction between these shells reads 

Uir3) = III ^D{r',)p'dV' . (4.1) 

Here, ra > i?2 + R2 is magnitude of three-dimensional vector between centers of shells 
and rJj is magnitude of three-dimensional vector between center of the first shell and 
an arbitrary elementary mass dm' = p'dV within second shell. Potential ^fDif'^) is 
given by Eq. (|3.7[) . After integration, we obtain: 

GNinm' GnPp' 
Uirs) = — — 

X V' ^e-2-X'^'-3[/i,.(i?2) - hk{Ri)][hk{R'2) - hk(R[)] . (4.2) 
k 

In the case of two spheres [Ri = 0,i?2 = R and R'^ = 0,i?2 ^ R'), this formula is 
reduced to: 

U{r,) = - ^ Y; \e^^^-^-h,mhk{R') . (4.3) 

rs 47r4r3 ^ Xk 
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For working approximation, it is often sufficient to keep in sums the first Kaluza-Klcin 
modes, i.e. to use the Yukawa approximation. In this approximation, (|4.2p and (|4.3p 
can be rewritten correspondingly: 

GNmm' IGn^GNPp'ct^^ -rs/x 



Uirs) 



and 



where 



h{R) 



GNfnm' 



1 + 9a ( - 



R' 



hml\h(R)Q^ (4.4) 
h{R)h{R'), (4.5) 



^cosh(^)-smh(| 



Here, A = a/(27r) and parameter a = 2s is described in p. 91 



(4.6) 



4-2. Gravitational force between two spheres 

In Yukawa approximation, gravitational force between two spheres is 
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dU 
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G]\rmm' 



l + 9a 



A \ / A \ rs 
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h{R) h{R') 



9 ^AV/^AV!^ p-('-3-«-«')A 



1 + -a , 
4 \R 



R' 



A 



(4.7) 



where in the last expression we use conditions R, R' A. If surfaces of the spheres 
are on the distances of the order of the maximal period: r^ — R — R' a, then we get 



^^(^3) 



GN'mm' 



1 + 0.0084s 



A 



Giqmm' 



(4.8) 



For example, in the case of one extra dimension and R = R' = 1cm, A = 4.7 x 10^'^cm, 
we get for deviation from the newtonian formula the estimate S = 1.8 x 10^*^. To obtain 
it, it is necessary to remember that for d = 1 parameter s = 1. 



5. Smeared extra dimensions 

Now, we get onto the asymmetrical extra dimension models (cf. with [5^ ) with topology 



D 



j,d~p ^j,p ^ p<d, 



(5.1) 



where we suppose that {d — p) tori have the same "large" period a and p tori have 
"small" equal periods b. In this case, the fundamental Planck scale relation (|2.5p reads 



(5.2) 



Additionally, we assume that test bodies are uniformly smeared/spreaded over small 
extra dimensions. Thus, test bodies have a finite thickness in small extra dimensions 
(thick branc approximation). For short, we shall call such small extra dimensions as 
"smeared" extra dimensions. If p = rf then all extra dimensions are smeared. 
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It is not difficult to show that the gravitational potential does not feel smeared 
extra dimensions. We can prove this statement by three different methods. First, 
we can directly solve D-dimensional Poisson equation (|2.ip with the periodic 
boundary conditions for the extra dimensions ^p4-i,...,^d and the mass density 
p = {m/ Yl^^i 6{r3)6{^p+i)...S{^d)- As a result, we obtain the unique solution 
V-D(r3, Cp+i; ■■■iS.d) of the form (j2.3p which does not depend on ^i, . . . , and satisfies 
the limit lim <^_D(r3, ^p+i, ^d) = 0. 

rs—t+oo 

Second, we can average solutions (j2.3p (where point mass m is replaced by 
mYli^iid^.i/ai)) over smeared dimensions ^i,...,'Cp and take into account that 
/g" cos(27rfc^/a)d^ = 0, a for fc 7^ and k = 0, respectively. Thus, the terms with 
these variables disappear from (|2.3p . For example, direct integration of the compact 
expression (j2.7p results in pure newtonian formula: 



smh 

Gnui 



cosh ( ^ I — cos I — ^ I d^ 



a 



arctan 



'cosh(2ip) + iy^' fn^ 



cosh(2zni) _ 1 



tan 



a 







^^"^ (5.3) 



Finally, from the symmetry of the model, it is clear that in the case of a test mass 
smeared over extra dimensions, the gravitational field vector E^j = —'^d^d does 
not have components with respect to smeared extra dimensions. For example, if we 
consider the model with topology of a cylinder S"^ x R where is the smeared extra 
dimension and R is the external dimension, then E2 will be parallel to an element of 
the cylinder. For simplicity, let us consider a model where all extra dimensions are 
smeared. Then, the Poisson equation (|2.ip can be rewritten in the form: 

Vd^d = -SdGvPd{yd) = -47rGArm(5(r3) , (5.4) 
where we use relation (12. 4p between gravitational constants and formula p£){v]j) ~ 

d 

m5{v^)l Y\ tti. After integrating both parts of this formula over multidimensional 
volume y = Va U Vintemai we obtain: 

/ a d°V = -4:nGNmV,nternal = COUSt . (5.5) 

Jv 

Here, d^V = £\ai,...,ao\d^"'^ A ... A da;"° — dxdydzdx^^ . . . dx'^'^ and ai = 
1, 2, 3, ?i, . . . , id- Indexes ii, . . . ,id enumerate extra dimensions. Applying the Gauss 
theorem, we get 

/ ad^'V^ f E'^^d^'-^J:^^ ^ AnRlEoVMernal , (5.6) 

JV ' ' JdV 

where 

d^-^S^,. = £a.|c,...oj rfx"^' A. . . Ada:"", (5.7) 

D - 1 

and in three-dimensional spherical coordinates d^^^E^a — R3 sin dd-ddipdx^^ . . . dx^°-^ . 
To get (|5.6p . we use the fact that i?^ = Ed is the only non-zero component of the 
gravitational field vector. Therefore, comparing (15. Sp and (|5.6p . we obtain: 



Eoirs) = — =^ VD{r^)^ . (5.8) 
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Thus, all these 3 approaches show that the gravitational potential does not feel 
smeared extra dimensions. It means that in the case of p smeared extra dimensions, 
the wave numbers fci , . . . , fcp disappear from (|2.3p and we should perform summation 
only with respect to fcp+i, . . . , fc^. 



5.1. Effective gravitational constant 

As it follows from (j4.7p . in the Yukawa approximation, the gravitational force between 
two spheres with masses mi, m2, radiuses i?2 and distance ra between the centers 
of the spheres reads: 

Fir,) = _^W^3)™ (3 



where 



' —\ !lp-('-3--R.i--R2)/A 



R2J 



= Gn{1 + Sg) (5.10) 

is an effective gravitational constant. Here, we took into account that in the case 
of p smeared dimensions, the prefactor a should be replaced by 2(d — p). Now, 
we want to evaluate the corrections 5q to the Newton's gravitational constant and 
to estimate their possible influence on the experimental data. As it follows from 
Figure 2 in the CODATA 2006, the most precise values oi Gn were obtained in 
the University Washington and the University Ziirich experiments [HI [TO] • They are 
GAr/lO-^m^kg^^s-^ = 6.674215 ± 0.000092, and 6.674252 ± 0.000124, respectively. 
Let us consider two particular examples: {T> = 5)-dimensional model with d — l,p = 
a — 2 and (T> — 10)-dimensional model with d — 6,p = S a — 6. For 
these values of a, Figure 6 in [7] gives the upper limits for A = a/(27r) correspondingly 
A « 4.7 X 10~'^cm and A ~ 3.4 x 10~^cm. To calculate Sq, we take parameters 
of the Moscow experiment [TT]: Ri w 0.087cm for a platinum ball with the mass 
mi — 59.25 X 10~'^g, R2 ~ 0.206cm for a tungsten ball with the mass m2 = 706 x lO^^g 
and r, = 0.3773cm. For both of these models we obtain Sq ~ 8.91 x lO^^'^ and 
Sq ~ 1.06 X 10~^^, respectively. These values are rather far from the measurement 
accuracy of Gn in [HI [TO]- In future, if we can achieve in the Moscow- type experiments 
the accuracy within these values, then, changing radii i?i^2 and distance r,, we can 
reveal extra dimensions or obtain experimental limitations on considered models. 
We should note that small changes in experimental bounds for A result in drastic 
changes of Sq. For example, if for A we take the upper limits A ~ 2 x 10~^cm and 
A « 1.3 X 10~^cm which follow from previous experiments (see Figire 5 in [I] for a = 2 
and 6, respectively) then Sq « 0.0006247 and So ~ 0.0000532 and these figures are 
very close to the measurement accuracy of Gat in [51 llOj. 



5.2. Model: P = 10 with d^6,p = 3 

Let us consider in more detail (P = 10)-dimensional model with 3 smeared dimensions. 
Here, we have very symmetrical with respect to a number of spacial dimensions 
structure: 3 our external dimensions, 3 large extra dimensions with periods a and 
3 small smeared extra dimensions with periods b. For b we put a limitation: 
b < bjnax = 10~^^cm which is usually taken for thick brane approximation. This 
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limitation follows from the electrical inverse-square law experiments. Although, in the 
next subsection we shaw that such approach can be significantly relaxed for models 
with smeared extra dimension, here we still use this bound. 

As we mentioned above, in the case of a = 6, for a we should take a limitation 
a < a.max = 2.1 X 10~^cm. To solve the hierarchy problem, the multidimensional 
Planck scale is usually considered from ITeV up to approximately 130 TeV (see e.g. 
[8l[12]). To make some estimates, we take Mmin = lTeV< Mp;(io) < M^ax = 50TeV. 
Thus, as it follows from (|5.2p . the allowed values of a and b should satisfy inequalities: 



5. MS ^ « " c MS . ■ ^ ^ ^ 

Counting all limitations, we find allowed region for a and b (shadow area in Figure 
HI). In this trapezium, the upper horizontal and right vertical lines are the decimal 
logarithms of amax and bmax, respectively. The right and left inclined lines correspond 
to Mp;(io) = ITeV and Mp;(io) = 50TeV, respectively. To illustrate this picture, 
we consider two points A and B on the line Mp;(io) — ITeV. Here, we have 
o = 2.1 X lO-^cm, b = 10-2°-9cm for A and a = IQ-^cm, b = IQ-^^-^cm for B. These 
values of large extra dimensions a are much bigger than in the standard approach 
a ~ 10(32/6)-i7cjn « 10-"- W nig. 



Ig a 




Ig b 



Figure 4. Allowed region (shadow area) for periods of large (a) and smeared 
(fe) dimensions in the model D = 10 with d = 6,p = 3. 



5.3. Model: V- arbitrary and d — p 

In this model, the test masses are smeared over all extra dimensions. Therefore, in non- 
relativistic limit, there is no deviation from the Newton's law at all. Surprisingly, this 
result does not depend on the size of smeared extra dimensions. The ISL experiments 
will not show any deviation from the Newton's law without regard to the size b (see 
also ()5.10p where d — p = 0). Similar reasoning are also applicable to Coulomb's law. 
It is necessary to suggest other experiments which can reveal the multidimcnsionality 
of our spacetime. Nevertheless, we can solve the hierarchy problem in this model 
because (j5.2p (where d = p) still works. For example, in the case of bosonic string 
dimension P = 26 we find A/p;(26) ~ 31TeV for b = 10~"'^^cm. In the case V = 10 
we get Mp;(io) « SOTeV for b = 5.59 x 10~^'*cm. In general, if we suppose that 
ITeV < A/pi(23) < aOTeV then we obtain: 



p=l, 1.65 X 10"cm < 6 < 4.46 X lO^^cm 
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p = 2, 3.81 X lO^'^cm < 6 < 3.43 X lO^^cm, 
p = 3, 4.83 X lO^^cm < 6 < 1.40 x lO^^cm, 

p = 6, 5.59 x lO^^Vn < 6 < 5.21 x lO^^^cm, 
where lower limit corresponds to 30TeV and upper limit agrees with ITeV. 

6. Conclusions 

We have considered generalization of the Newton's potential to the case of extra 
dimensions where multidimensional space has topology Md = M.'^ x T"^ . We obtained 
the exact solution (j2.3p which describes smooth transition from the newtonian 
behavior 1 /ra for distances bigger than periods of tori (the extra dimension sizes) to 
multidimensional behavior 1 /^g^^ in opposite limit. In the case of one extra dimension, 
the gravitational potential is expressed via compact and elegant formula (|2.7p . We 
applied these exact solutions to some practical problems to get the gravitational 
potentials for considered configurations. Found potentials were used to calculate the 
acceleration for point masses and gravitational self-energy. 

To estimate corrections to the newtonian expressions, it is sufficient to keep 
only the first Kaluza-Klein modes. Then, we obtained that the correction term 
has the form of the Yukawa potential with parameters defined by multidimensional 
models. Such representation gave us a possibility to use the results of the inverse 
square-law experiments to get limitations on periods of tori. As an Yukawa potential 
approximation, it was shown that in the Cavendish-type experiments the corrections 
(due to the extra dimensions) to the newtonian gravitational constant are still far 
from the measurement accuracy of experiments for the determination of Gat. 

Then, we proposed models where test masses can be smeared over the extra 
dimensions. The number of the smeared dimensions can be equal or less than the 
total number of the extra dimensions. We proved that the gravitational potential 
does not feel the smeared dimensions and this conclusion does not depend on the size 
of these dimensions. Such approach opens new remarkable possibilities. For example 
in the case P = 10 with 3 large and 3 smeared extra dimensions and Mp;(io) = ITeV, 
the large extra dimensions can be as big as the upper limit established by the ISL 
experiments for a — 6, i.e. a « 2.1 x 10~^cm. This value of a is in many orders 
of magnitude bigger than the rough estimate a « 10^^^'^cm obtained from the 
fundamental Planck scale relation of the usual form ()2.5|) . The limiting case where 
all extra dimensions are smeared is another interesting example. Here, there is no 
deviation from the Newton's law at all. Thus, these models can explain negative 
results for the detection of the extra dimensions in ISL experiments irrespective of the 
size of the extra dimensions. Nevertheless, these models can be still used to solve the 
hierarchy problem. 
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Appendix A. The gravitational force law under dimensional reduction 

Let us consider a product space AId consisting of two components Md = M3 x 

Md, where for simplicity of the subsequent calculations we assume M3 = IR'^ 

as 3— dimensional flat external space and a d— torus Md — T'^ (with the same 

characteristic length a along each of the d dimensions) as compact internal space. 

The volume of Md is hence given as vol{Md) = Vd = a'^. 

The aim of this appendix is to demonstrate that the force laws for the full theory 

on AI]j and the effective theory on M3 have correspondingly the form 

p. . ^ TO1TO2 
F{r3+d) = Gt> ^ 

Fir3) = G,^, (A.2) 
^3 

where gravitational constants Gt> and Gn are related with each other in accordance 
with (12. 4p . To perform it, we consider separately the regimes of small distances 
''a+d ^ compared to the size of the internal space and of large distances 3> a. For 
f3+d ^ a the field lines penetrate all D space dimensions, whereas at length scales 
ra ^ a the internal dimensions are not accessible and effectively the field lines are 
spreading only along the 3 spatial dimensions of the external space M3. 

As in [3], we pass for simplicity of the calculations from the compact d— torus to 
its covering space M''. The single mass mi in T"^ is then mapped to a mass lattice in 
IR'*, which is build from rrii and its mirror images in the cover. 

For a test mass TO2 at a small distances rs+d <C a from mi the mirror masses give 
a negligible contribution to the gravitational force and in first order approximation 
the Gauss law in D space dimensions holds in accordance with (jA.l[) . In this 
case the infinite set of mirror masses can in rough approximation be considered as 
symmetrically distributed around mi, m2 so that their contributions will compensate 
each other. 

For the analysis of the opposite case of a large separation ^ a between the 
masses mi and m2, we assume the mass mi placed in the origin of the Z?— dimensional 
space. The lattice of the mirror masses spans then a d— dimensional subspace M'' 
("wire") in the space M"* x M'^ and the total effective force F^ff will result from the 
gravitational attraction of the whole lattice-subspace. For rs ^ a, the mirror masses 
appear homogeneously smeared about this subspace with a "surface" density rrii/a'^. 
Let us choose coordinates 

(z,r3) e M'' X (A.3) 

and a test mass m2 located at the point (0,r3). The force orthogonal to the lattice 
produced on m2 by a small volume element dVd at a point (z, 0) of the (lattice) covering 
space M'' is given by 

r3 midVd 1712 



dF{z,r3) Gv^^^^^2y/2 ad (^2+^2)(2+d)/2 

Gvmim2 rsdVd 

(z2 ^^2y3+<i)/2 • 



(A.4) 
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The factor r-iji^z^ + fg)^/^ is the cosine between the force direction and the lattice 
normal. Due to the lattice symmetry, the force component parallel to the lattice 
is compensated by an opposite force component from the point (— z, 0). The total 
effective force ^"e//(?'3) of the mass lattice on the test particle can now be easily 
obtained by integrating over the volume of M'^. Choosing spherical coordinates on R'', 
the volume of a thin shell with radius z is 

dVi = Siz'^-^dz, (A.5) 

where Sd is square of [d — l)-dimensional sphere of a unit radius (see (|2.1[) ) and the 
total effective force can be calculated as 

Gvmxvni f°° z'^-^dz 
With the substitution t ^ z^ the integral can be transformed to the standard integral 

poo j-a-lj^ 

/ — =b°'~PB(a,p-a), 0<5ftQ<5ftp (A.7) 

Ja + by 

so that for the corresponding term in (|A.6p we find 

z^~^dz 1 r°° f^i^-^t 



(z2 + r2)(3+'i)/2 2 7o (t + r2)(3+'i)/2 

_ I „(di\_ i r(f)r(|) 

Expressing the F— functions in (|A.8P in terms of surface areas: r{D/2) — 2tt^/^/Sd 
gives 



-irfz 1 S- 



^^'^ (A.9) 



(z2 + r2)(3+d)/2 rlSdSa' 
Hence, the total effective force takes the simple form 

^^e//(r3)^^^^ (A.IO) 

and equating it with the gravitational force law (jA.2|l in three-dimensional space, 
we reproduce the result (j2.4|) for the relation between three-dimensional and 
multidimensional gravitational constants. 
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